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The hot and dense matter generated in heavy-ion collisions may contain domains which 
are not invariant under P and CP transformations. Moreover, heavy-ion collisions can gen¬ 
erate extremely strong magnetic fields as well as electric fields. The inferplay between the 
electromagnetic field and friangle anomaly leads fo a number of macroscopic quanfum phe¬ 
nomena in these P- and CP-odd domains known as the anomalous transports. The purpose 
of the article is to give a pedagogical review of various properties of the electromagnetic 
fields, the anomalous transports phenomena, and their experimental signatures in heavy- 
ion collisions. 
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I. INTRODUCTION 

As is well known, the strong interaction provides the mechanism that binds the quarks and 
gluons together to form the hadrons such as the proton and neutron. Our contemporary under¬ 
standing of strong interaction is described by quantum chromodynamics (QCD) — a quantum 
gauge field theory based on color S Lf (3) gauge symmetry. Despite its simple form, QCD possesses 
a number of remarkable properties among which the most mysterious one may be the color con¬ 
finement: in vacuum, the quarks and gluons, as being colorful particles, are always confined in 
colorless hadrons. The color confinement property forbids us to observe isolated quark or gluon. 
However, at high temperature and/or high quark chemical potentials, the normal hadronic mat¬ 
ter is expected to transform to deconfined quark-gluon matter. When the temperature is high 
such quark-gluon matter is usually called the quark-gluon plasma (QGP). Lattice QCD simula¬ 
tions which are the first-principle computations that solve QCD directly show that the "transition 
temperature" from hadronic matter to QGP is about 200 MeV (at zero quark chemical potentials). 
It is believed that such a high temperature was once realized in the universe at a few microsec¬ 
ond after the Big Bang and the QGP was created at that time. In laboratory, up to now, the only 
method to achieve such a high temperature is to use the high-energy heavy-ion collisions. Such 
collisions have been carried out in the Relativistic Heavy Ion Gollider (RHIC) at Brookhaven Na¬ 
tional Laboratory (BNL) and in the Large Hadron Gollider (LHG) at the European Organization 
for Nuclear Research (GERN). The top center-of-mass energy per nucleon pair at RHIG Au + Au 
collisions is \/s = 200 GeV and at LHC Pb + Pb collisions is \/s = 2.76 TeV and will be upgraded 
to \/s = 5.5 TeV soon. In these colliders, two nuclei are accelerated to velocity very close to the 
speed of light and then collide, the energy deposition in the reaction region can be large enough 
to create the hot and dense environment in which the deconfinement condition is reached. Mea¬ 
surements performed at RHIC and LHC have collected many signals supporting the generation 
of the QGP and have also revealed a variety of unusual properties of the QGP, e.g., its very small 
shear viscosity comparing to the entropy density and its high opacity for energetic jets. 

Heavy-ion collisions can generate electromagnetic (EM) fields as well [1]. Recent numeri¬ 
cal simulations found that the magnitude of the magnetic field in RHIG Au + Au collisions at 
y/s = 200 GeV can be at the order of 10^® — 10^^ Gauss ^ and in LHG Pb + Pb collisions at 

^ In the relevant literature, people usually use or MeV^ as the unit of eB where e is the electron charge magnitude 
and m„ fa 140 MeV is the pion mass. In converting to the SI or Gaussian units, it is helpful to note the following 
relation: 1 MeV^ = e • 1.6904 x 10^^ Gauss (if = c = 1, otherwise the right-hand side should be multiplied by hc^). 
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^/s = 2.76 TeV can reach the order of 10^*^ Gauss [2-11], The electric filed can also be gener¬ 
ated owing to event-by-event fluctuations [5, 7-9] (we will explain the physical meaning of such 
fluctuations in Sec. II) or in asymmetric collisions like Cu + Au collision [12-14], and its strength 
is roughly of the same order as the magnetic field. Thus heavy-ion collisions provide a unique 
terrestrial environment to study QCD matter in strong EM fields. In particular, recently it was pro¬ 
posed that the magnetic field can convert topological fluctuations in the QCD vacuum into global 
electric charge separation along the direction of the magnetic field. The underlying mechanism is 
the so-called chiral magnetic effect (CME) [2,15]. Some relatives of the CME were also proposed, 
including the chiral separation effect (CSE) [16, 17], chiral electric separation effect (CESE) [18], 
chiral magnetic waves (CMW) [19], chiral vortical effect (CVE) [20-22], chiral vortical wave [23], 
chiral heat wave [24], chiral Alfven wave [25], etc. They all represent special transport phenom¬ 
ena that are closely related to chiral anomaly ^ and thus are called anomalous transports. The 
experimental searches of the anomalous transports have been carried out at RHIC and LHC and 
the measurements indeed offered signals consistent with the predictions of the CME, CMW, and 
CVE; see the discussions in Sec. IV. 

The purpose of this paper is to give a pedagogical review of recent progresses on the study of 
EM fields and the anomalous transport phenomena induced by EM fields in heavy-ion collisions. 
We will keep all the discussions as intuitive as possible and lead the readers who wish to under¬ 
stand more technical details to proper literature. To access this paper, the readers do not need to 
have expert knowledge of QCD; only elementary knowledge of quantum field theory and heavy- 
ion collisions are needed. (Perhaps the only exception is Sec. Ill A (3) where some knowledge of 
topology and gauge field theory is needed; we thus give more thorough discussion and put nec¬ 
essary references there so that the readers can easily trace the relevant literature.) Thus this paper 
will be particularly useful for graduate students who have finished their first-year courses and 
wish to enter the exciting research area of anomalous transport phenomena. In this aspect, the 
present review is complementary to existing excellent reviews, e.g.. Refs. [26-31] in which more 
advanced materials can be found. 

We organize the paper as follows. In Sec. II, we shall discuss some general properties of the 
EM fields in heavy-ion collisions. In Sec. Ill, we shall give an elementary introduction to the 
anomalous transports in parity-odd (P-odd) and/or charge-conjugation-odd (C-odd) medium. 
The experimental implications of the anomalous transports and the current status of their detec¬ 
tion in heavy-ion collisions will be reviewed in Sec. IV. Some discussions will be presented in 
Sec. V. 

In addition to the anomalous transports, the EM fields can drive a range of other intriguing 
phenomena including, for example, the magnetic catalysis of chiral symmetry breaking [32-34], 
the inverse magnetic catalysis or magnetic inhibition at finite temperature and density [35M9], the 

^ We will not distinguish the terms "chiral anomaly", "axial anomaly", and "triangle anomaly" in this article. 
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possible p meson condensation in strong magnetic field [50-54], the neutral pion condensation in 
vaccum [55], the anisotropic viscosities in hydrodynamic equations [56-60], and the early-stage 
phenomena in heavy-ion collisions like the EM-field induced particle production [26, 61-66] and 
the dissociation of heavy-flavor mesons [67-71], These topics will not be the main focus of this 
article. Some of them are nicely reviewed in Refs. [28, 72, 73]. 


II. PROPERTIES OE ELECTROMAGNETIC FIELDS IN HEAVY-ION COLLISIONS 

The reason why heavy-ion collisions can generate magnetic fields is simple: nuclei are pos¬ 
itively charged and when they move they generate electric currents which in turn induce the 
magnetic fields. In a noncentral heavy-ion collision, two counter-propagating nuclei collide at a 
finite impact parameter b; one can easily imagine that the magnetic field at the center of the over¬ 
lapping region will be perpendicular to the reaction plane owing to the left-right symmetry of the 
collision geometry (see Fig. 1 for illustration). However, in a real collision event, this left-right 
symmetry may be lost because the nucleon distribution of one nucleus would not be identical to 
another. We will come to this point later, but first let us estimate how strong the magnetic field 
can be. 

Let us consider Au + Au collisions at fixed impact parameter b = 10 fm and at RHIC energy 
^/s = 200 GeV as an example. If we approximate the problem by assuming that all the protons 
are located at the center of the nucleus, then by naively applying the Biot-Savart law we obtain 

- eBy ~ 2Zau7^?^z Q j ~ ~ 10^^ Gauss, (1) 

where Vz = vl —(2m]vVA^s)2 ~ 0.99995 is the nucleon mass) is the velocity of the nucleus, 
7 = Ij ^J\ — v\ ~ 100 is the Lorentz gamma factor, and Zau = 79 is the charge number of gold 
nucleus. The minus sign on the left-hand side is because the magnetic field in pointing to the —y 
direction in the setup shown in Fig. 1. 

This is really a huge magnetic field. It is much larger than the masses squared of electron, 
and light quarks, and thus is capable of inducing significant quantum effects. It is 

also larger than the magnetic fields of neutron stars including the magnetars which may have 
surface magnetic fields of the order of 10^^ — 10^^ Gauss [74, 75]. Therefore the magnetic fields 
generated in high-energy heavy-ion collisions are among the strongest ones that we have ever 
known in current universe. (In the early universe, there was a possibility to generate an even 
stronger magnetic field through the electro-weak transition, see Ref. [76] for review.) One can 
expect that such a huge magnetic field may have important consequences on the dynamics of the 
quark-gluon matter produced in heavy-ion collisions. We will discuss several such consequences, 
namely, the anomalous transport phenomena in Sec. III. In this section we will focus on the fields 
themselves. 
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FIG. 1. The geometrical illustration of the noncentral heavy-ion collision. ITere b is the impact parameter 
and is the radius of the nucleus. The magnetic field B is expected to be perpendicular to the reaction 
plane due to the left-right symmetry of the collision geometry. 

A. Computations of the electromagnetic fields 


The estimation given in Eq. (1) is too simplified, in order to have a more reliable simulation 
for the electromagnetic (EM) fields in heavy-ion collisions, the following issues need to be taken 
into account. (1) We need the knowledge of the proton and neutron distributions in a given 
nucleus. Eor this purpose, we can choose the well established Woods-Saxon distribution to use. 
(2) In real heavy-ion collisions, because the proton distribution varies from one nucleus to another, 
the generated EM fields vary from event to event. It is thus important to study the event-by- 
event fluctuation of the EM fields [5, 7-9]. (3) We need to replace the Biot-Savart law by the full 
relativistic Lienard-Wiechert potentials which contain the retardation effect. 


eE{t,r) = 
eB{t,r) = 

n 


Rn RfiVii , 2 \ 

{Rn-Rn-Vny^ 


Vn y. Rn 
{Rn Rn ' Vny 


(1 


( 2 ) 

(3) 


where the summation is over all the charged particles, Z„ is the charge number of the nth particle, 
R„ = r — r„ is the relative position of the field point r to the source point r„ of the nth particle, 
Vn is the velocity of nth particle at the retarded time t„ = t — \r — r„\. Note that Eqs. (2)-(3) have 
singularities at Rn = 0; in practical calculations a variety of regularization schemes have been 
used and consistent results are obtained after taking the event average [3-5, 7-9, 13]. 

As the EM fields in heavy-ion collisions can be much larger than the electron and light quark 
masses squared, one may worry about the possible quantum electrodynamics (QED) correction to 
the otherwise classical Maxwell field equations (the Lienard-Wiechert potentials are the solutions 
of the Maxwell equations). So let us make a magnitude estimate of such QED correction by using 
the one-loop Euler-Heisenberg effective lagrangian for soft photons (see Ref. [77] for review): 




ds 
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where /f' is the electric current, A^‘ is the EM potential, is the strength tensor, F^^ = 

and X = \F^yF^^ + ^F^yF^-^. At strong-field limit, the asymptotic form of Eq. (4) behaves 

like [77] 


>Ceh ~ —^iiF 


1 - 


247r2 


In 


^|f2| 


m-i 


F^VF 


]lVr 


( 5 ) 


where F^ = F'^^F^^. The field equations derived from this lagrangian can be regarded as the 
Maxwell equations but with a renormalized charge (keeping only the leading-log term because 
|eF| > ml) 


e —)■ e 


1 - 


24712 


In- 


^|f2| 


-1 -1 


m2 


( 6 ) 


at leading-log order. Thus, we can find that even for very strong EM field, e.g., |eF| 100m2, the 
quantum correction can only amend the final restuls by a few percent. This justifies the applica¬ 
bility of Eqs. (2)-(3). 

In the following subsections we will review the recent results of the EM fields in heavy-ion 
collisions obtained by using the Lienard-Wiechert potentials on event-by-event basis. We will 
mainly focus on Au + Au collisions at RHIC and Pb + Pb collisions at LHC; other collision systems 
will be briefly discussed in Sec. IIH. 


B. Impact parameter dependence 

We first show in Eig. 2 the impact parameter dependence of the EM fields at r = 0 and 1 = 0 
where the initial time 1 = 0 is set to be the time when the two colliding nuclei completely overlap. 
The curves with full dots are for Au -t- Au collision at RHIC energy ^/s = 200 GeV and the curves 
with open dots are for fields scaled by a factor 13.8 = ^Slhc / V^RHIC for Tb + Pb collision at 
LHC energy -^/s = 2.76 TeV. In these figures (and also in the figures hereafter) (• • ■) represents the 
average over events. One can find that: 

(1) The event averaged EM fields have only one nonzero component, (By) ^ 0; all other compo¬ 
nents vanish: (B^) = (Bz) = (E) = 0. However, owing to the fluctuations of the positions of 
protons in the nuclei, their magnitudes in each event can be large (except for the z-components 
which are always small). This is reflected in the averaged absolute values of the fields and is most 
evident for central collisions [5, 7]. 

(2) When b < 2 Ra with the nucleus radius, the event-averaged field e{By) is proportional to b 
and it reaches its maximum value around IR^ ■ The fluctuation-induced fields are not sensitive to 
b when b < 2 Ra- 
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b(fm) 


FIG. 2. The EM fields at f = 0 and r = 0 as functions of the impact parameter b. Figures are modified from 
Ref. [7]. 


C. Collision energy dependence 


As investigated in Ref. [5, 7], to high precision, the magnitudes of EM fields linearly depend 
on the collision energy ^/s. Actually the absolute values of the EM fields satisfy very well the 
following scaling law, e|Field| oc ^/sf{b/R/^) where f{b/RA) is a universal function which has 
the shapes for {e\Bxy\) and (e|£x,i/|) as shown in Fig. 2. For the event-averaged magnetic field, 
e{By), the following formula approximately expresses its impact parameter b, collision energy 
y/s, charge number Z and atomic number A dependence: 


e{By) 


y/s Z b 2 

2m^ A2/3 2R,4™” 


for b < 2 Ra- 


( 7 ) 


Note that the prefactor y/s/ {2mjy) is nothing but the Lorentz gamma factor. 


D. Spatial distributions 

The spatial distributions of the EM fields are evidently inhomogeneous. The contour plots for 
{eBx,y,^, {eEx,y^z)f {^\Bx,y,z\)f and (e|Ex,y,z|) in the transverse plane at impact parameter & = 10 fm 
and at f = 0 for RHIC energy are shown in Fig. 3. The distributions of the fields for LHC energy 
is merely the same but with 2760/200 = 13.8 times larger magnitudes according to Sec. IIC. 

One may notice that for noncentral collision, the y-component of the electric field is very large 
along the y-direction, reflecting the fact that at f = 0 a large amount of net charges stays tempo¬ 
rally in the center of the "almond"-shaped overlapping region. This strong, out-of-plane electric 
field may drive positive (negative) charges to move outward (toward) the reaction plane, and 
thus induce an electric quadrupole moment in the produced quark-gluon matter. Such an electric 
quadrupole moment, as argued in Ref. [78], may lead to an elliptic flow imbalance between 7r+ 
and TT^; see Refs. [7, 79] for the detail. 
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FIG. 3. (Color online) The spatial distributions of the EM fields in the transverse plane at t 
at RHIC energy. The unit is m^. The dashed circles indicate the two colliding nuclei. 
Ref. [7]. 


= 0 for b = lOfm 
Figures are from 


E. Azimuthal correlation with the participant planes 


We have seen that the event-by-event fluctuations of the nuclear distribution can strongly mod¬ 
ify the magnitudes of the EM fields, one then may ask how these event-by-event fluctuations af¬ 
fect the azimuthal orientations of the EM fields (see the illustrating Eig. 4). In fact, as revealed 
recently [8], the event-by-event fluctuations generally make the magnetic fields unaligned with 
the normal direction of the reaction plane. As a consequence of the event-by-event fluctuations 
of nuclear distribution, the distribution of the participants (the nucleons that participate in the 
collision) in the overlapping region varies from event to event as well. Thus for each event, the 
overlapping region is not perfectly almond-shaped and its short-axis may be rendered away from 
the impact-parameter direction. Such shape and direction variations can be captured by the so- 
called eccentricity parameters €„ and harmonic angles n = 1,2,3, . Mathematically, they 

are defined as 


lY J (frp{r)r^e“t’ 

€ie ^ = — 


= - 


f d^rp(r)r^ 


, n > 1, 


( 8 ) 

( 9 ) 


f d^rp(r)r’^ 

where p(r) is the transverse distribution function of the participants. If there were no event-by- 
event fluctuations, Y 2 for each event should be equal to Y^p. As we will see, the azimuthal ori¬ 
entation between ips of B (and tpE of E) and Y 2 (the second harmonic angle of the participants) 
fluctuates with sizable spread in their relative angle tpB — Y 2 about the expected value nil.. (Note 
that this would imply that i/’s — Yrp and i/’s — Yp with Yrp and Yp the reaction plane angle and 
event plane angle also fluctuate.) This can be clearly seen from Eig. 5 which shows the histograms 
of — Y 2 over events for different h. Eor the most central case the events are uniformly dis¬ 
tributed indicating negligible correlation between ips and Y 2 ; while for noncentral collisions the 
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event distribution behaves like a Gaussian peaking at xps —^2 — ti/T. indicating correlation be¬ 
tween Yb and Y 2 . For the following reason, such fluctuations in the correlation between tps 
(as well as ips) and the participant planes may have important impacts on the experimentally 
measured quantities. 



FIG. 4. Illustration of the azimuthal angles of the participant plane Ypp and electric and magnetic fields. 


Let us consider the chiral magnetic effect (CME, see Sec. Ill A) as a concrete example, but 
the analysis can be extended to the observables of other EM-field-induced transport phenomena. 
The CME contributes to the single particle distribution for charged hadrons a component oc 
qeB cos{(p — ip b) with q = Y the charge which in turn contributes to the two-particle distribution 
the following term, 

fcjoccip 0 ^ qccq[i{eBfcos{(poc - iPb) cos (^^ - ipB), ( 10 ) 

where (pa^^ are the azimuthal angles of the hadrons a and /3. We therefore can extract the CME 
contribution to the two-particle correlation 7 ^^ = {cos{(poc Y(p^~ 2 ^ 2 )) (which was used by STAR 
and ALICE collaborations to detect the charge separation with respect to the reaction plane, see 
Sec. IV A) as 

locfi qccqfi{{eBfcos[2{ipB-Y2)]). (11) 

If the S-direction were always perpendicular to the reaction plane while Y 2 always coincide with 
Yrp (which we set to be zero here), then we simply have 7 a^ oc —qaq^{{eB)^). But the fluctuations 
in magnetic field as well as in participant planes will blur the relative angle between the two and 
modify the signal by a factor ~ {cos[2{ipB — Y 2 )]). (Here we note that the magnitude of the mag¬ 
netic field has no noticeable correlation to its azimuthal direction [ 8 ].) Similarly, if one measures 
the two-particle correlation with respect to higher harmonic participant plane, for example, the 
fourth harmonic plane Y 4 , (cos[2(^a. + (p^ — 2 Y 4 )]), the azimuthal fluctuations of B will again 
contribute a modification factor ~ (cos[ 4 (i/ 7 b — Y 4 )]) to it. 

In Fig. 6 (left panel) one can find the computed average values of {cos[n{ipB — Y„)]) as func¬ 
tions of the impact parameter from the event-by-event determination of the S-field direction xp b 
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FIG. 5. The event-by-event histograms of — Y 2 at impact parameters b = 0,5,10 ,12 fm for Au + Au 
collision at RHIC energy. Here is the azimuthal direction of B field (at f = 0 and r = (0,0,0)) and Y 2 
is the second harmonic participant plane. This figure is from Ref. [8]. 



FIG. 6. The correlations {cos[n{xpB — Y,,)]) and {cos[n{xpE — Y„)]) for B- and £J-fields at the center of the 
overlapping region, r = 0, as fimctions of impacf paramefer. Figures are modified from Ref. [8]. 


at the collision center r = 0 and the participants harmonics, Y„, n = 1,2,3,4. The plots suggest: 

(1) The correlations between tps and the odd harmonics Yx,Y 3 are practically zero (as a conse¬ 
quence of parity invariance), while the correlations ofxpB with even harmonics Y 2 , Y 4 are nonzero 
but suppressed comparing to the non-fluctuating case. 

( 2 ) The centrality dependence of (cos[2(i/?b — Y 2 )]) agrees with the patterns shown in the his¬ 
tograms Fig. 5: it is significantly suppressed in the most central and most peripheral collisions 
indicating no correlations between xps and Y 2 while is still sizable for moderate values of b. 

(3) As checked in Ref. [ 8 ], there is no visible difference between the (eS)^-weighted correla¬ 
tion {{eB)^ cos[n{ipB — Y„)])/((eS)^) and the unweighted correlation {cos[n{ipB — Y,,)]) for 
n = 1,2,3,4. This indicates that the magnitude of the magnetic field does not noticeably cor¬ 
relate to its azimuthal direction. 

The event-by-event fluctuations also bring modification to the correlations between Fi-field 
orientation and the participant planes. In parallel to the S-field case, in Fig. 6 (right panel), we 
show the correlations {cos[n{ipE — Y„)]), n = 1,2,3,4, as functions of the impact parameter. It is 
seen that: 

(1) There is a sizable negative correlation (i.e., back-to-back) between xpE and Yi which is strongest 
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in the central collisions. This is simply because the pole of Yi with more matter will concurrently 
have more positive charges from protons which induce the £^-field pointing opposite to Yi. 

(2) There is also a weak correlation between tpE and Y 3 . 

(3) Similar to the S-field case, (e£^)^-weighted correlations {{eE)^cos[n{xj)E ~ '^n)])/{{sE)^) 
have no visible difference from the unweighted correlations {cos[n{xpE — Y„)]) indicating no cor¬ 
relation between £'-field magnitude and orientation. 

F. Early-stage time evolution 



FIG. 7. The early-stage time evolution of the electromagnetic fields at r = 0 with impact parameter b = 10 
for Au + Au collision at -y/s = 200 GeV and Pb + Pb collision at -y/s = 2.76 TeV. For t < lOfg, Eq. (13) fit 
the curve for (eBy) very well. After that, the remnants essentially slow down the decays of the transverse 
fields and Eq. (13) does not work well any more. Figures are modified from Ref. [7]. 


In a high-energy heavy-ion collision, right after the collision, the produced partonic matter is 
mainly consist of gluons and is in a far-from-equilibrium state. This partonic matter subsequently 
evolves toward thermal equilibrium and a large number of quarks and anti-quarks are excited 
during this thermalization process. Although so far we still lack a theory to quantitatively un¬ 
derstand the thermalization problem, the phenomenological studies revealed that the time scale 
of the completion of the thermalization is very short comparing to the total lifetime of the ther- 
malized quark-gluon plasma (QGP). The relevant information can be found in the review articles. 
Refs. [80-83]. Once the thermalization is locally achieved, the bulk evolution of the system can 
be well described by hydrodynamics. One of the transport coefficients of the hydrodynamics, 
namely, the electric conductivity cr has been numerically simulated by using lattice QCD recently 
and it was found that a for QGP is very large (see next subsection). A large a makes the QGP 
sensitive to the variation of the EM fields and which in turn strongly influence the time evolution 
of the EM fields themselves. Thus the time evolution of the EM fields in the QGP stage need spe¬ 
cial treatments and we leave this issue to next subsection. In this subsection we will focus on the 
stage before the thermalization is achieved (we call this stage the "early stage"). The quark-gluon 
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matter in the early stage is expected to be much less conducting than that in the QGP stage and 
we will simply assume it is insulating and thus ignore the response of the matter to the EM fields. 

In Fig. 7 one can find the numerical results of the early-stage time evolutions of the EM fields at 
r = 0 in collisions with b = 10 fm for Au + Au collision at = 200 GeV and for Pb + Pb collision 
at y/s = 2.76 TeV. The results are from Ref. [7]. The contributions to the EM fields come from the 
charged particles in spectators, participants, and remnants. We can see that the transverse fields 
decay very fast after the collision reflecting the fact that the spectators are leaving the collision 
region very fast. Once the spectators are all far away from the collision region, the remnants 
which moves much slower than the spectators become important and they essentially slow down 
the decays of the transverse fields. The lifetime of the magnetic fields due to spectators can be 
estimated as 


ts ^ Ra/{7Vz) ^—^Ra, ( 12 ) 

which is just half the time that one proton needs to pass through the nucleus freely. The lifetime 
fg is very short for large y/s: for Au + Au collision at 200 GeV, fg ~ 0.065 fm, while for Pb + Pb 
collision at 2.76 TeV, fg ~ 0.005 fm. Within the time period (0, fg) the fields decay slowly and 
when t > fg, we can approximate the early-stage time evolution of the event-averaged magnetic 
field, that is, {eBy) as 




{eByjO)) 

(l + f2/f|)3/2- 


(13) 


This formula works better for larger impact parameter b and larger y/s. As seen from Fig. 7, 
Eq. (13) fit the simulation results for {eBy{t)) very well for time f < lOfg; after that the remnants 
dominate and the separation between the curves from Eq. (13) and from the simulations are visi¬ 
ble. The Eq. (13) shows that for f > fg the magnetic field decays fast, {eBy{t)) ~ fg/1^ However, 
if at that time fg the QGP has been already formed, its EM response will significantly modify the 
time evolution of the fields. 


G. Late-stage (QGP-stage) time evolution 

The discussions and simulations presented in the last subsection are based on the assumption 
that the produced matter is ideally insulting. This assumption is adoptable only in the early stage 
where the system is gluon-dominated but becomes less and less justified as the system evolves 
and more and more quarks and anti-quarks emerge. As a matter of fact, the QGP is a good conduc¬ 
tor according to the theoretical and lattice QGD studies. At very high temperature the perturba¬ 
tive study gives that the electric conductivity of QGP is c « 6T/e^ [84]. An old lattice calculation 

^ One should note that the magnitude of the magnetic field is still very large even at 1 > lOlg; for example, (eBy) ~ 40 
MeV^ for RHIC Au + Au collisions at 1 = 1 fm which is still comparable to the light quark mass squared, 
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with Nf = 0 found that u ~ 7 Cem^ [85] at 1.5Tc < T < 3Tc with Tc the deconfinement tempera¬ 
ture. Another quenched lattice simulation using staggered fermions found that a « 0.4Cem7’ [86 ]. 
Recent quenched lattice studies using Wilson fermions obtained that a « (l/3)CEMr — CemI" [87- 
89] for temperature Tc < T < 3Tc. The lattice calculation with Nf = 2 dynamical Wilson fermions 
found that a ^ 0.4CemT at T ~ 250 MeV [90]. Another new lattice simulation using Nf = 2 T 1 
fermions obtained that cr ~ O.ICemT — 0.3CemT for temperature Tc < T < 2Tc [91, 92]. In these 
results, the EM vertex parameter Cem = YL/^ff and cjf is the charge of quark with flavor /; for 
example, Cem = (5/9)e^ if u,d quarks are considered while Cem = (2/3)e^ if u,d,s quarks are 
considered. Note that the deconfinement temperature Tc is different in = 0 and Nf ^ 0 cases; 
for example, if Nf = 0 we have Tc ~ 270 MeV while if = 2 -|- 1 we have Tc ~ 170 MeV. 

At T = 350 MeV and choosing a « 0.3CemT with u, d, s quarks contributing to Cem/ one 
can find that the resulted cr is about 10 ^ times larger than that of copper at room temperature 
(u-cu ~ 4.43 X 10-3 MeV at T = 20‘’C). 

Now let us analyze how the large cr influences the time evolution of the EM fields in the QGP 
stage which we refer to as "late stage". 

Our discussion will be based on magnetohydrodynamics. We first write down the Maxwell's 
equations, 

^ R 

VxE = -—, (14) 

dE 

VxB = — + J, (15) 

where J is the electric current. We treat the QGP as being locally charge neutral but conducting, 
thus J is the sum of the external one and the one determined by the Ohm's law, 

J = cr (E/ -|- u X B) Jextr (16) 


where v is the flow velocity of QGP and Jext is the current due to the motion of unwounded 
protons (most are in spectators). Using Eq. (16), we can rewrite the Maxwell's equations as mag¬ 
netohydrodynamic equations 


^ = Vx(uxB) + - (v^B 
dt a \ 

dE dv ^ 1 

— + — xB = vx{VxE) + - 
dt dt cr 


d^B 

'W 

V^E 


T V X dextj / 
Jext\ 

dfi dt ) ' 


(17) 

(18) 


where we have used the Gauss laws V ■ B = d and N ■ E = p = 0. Equation (17) is the induction 
equation, which plays a central role in describing the dynamo mechanism of stellar magnetic field 
generation. The first terms on the right-hand sides of Eqs. (17)-(18) are the convection terms, while 
the remained terms are called "diffusion terms" although they are not exactly in the diffusion- 
equation type. Let us discuss some outcomes of these magnetohydrodynamic equations. 
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(1) If = 0, that is, if the QGP does not flow, the Eq. (17) reduces to 


dB 

dt 


- - 


d^B 

'W 


-|- V X i7j 


ext 


(19) 


This equation can be solved by using the method of Green's function, the details can be found in 
Refs. [26, 93-97] in which the authors studied how the spectators induced magnetic field evolve 
in QGP phase (assuming the system is already in the QGP phase at the initial time). The main 
information from these studies are that the presence of the conducting matter can significantly 
delay the decay of the magnetic field. This is easily understood as the consequence of the Faraday 
induction: a fast decaying external magnetic field induces a circular electric current in the medium 
which in turn causes a magnetic field that compensates the decaying external magnetic field. 

For late times, the external current Jext from fhe spectators can be neglected (we will always 
assume this case in this subsection hereafter). If u 3> 1 /fc wifh tc the characteristic time scale over 
which the field strongly varies, one can neglect the second-order time derivative term and render 
Eq. (19) a diffusion equation: 


dB 

dt 


-V^B. 

a 


( 20 ) 


This case was studied in Ref. [61]. This equation describes the decay of the field due to diffusion, 
and the diffusion time of the magnetic field is given by 


to = L^cr, 


( 21 ) 


with L a characteristic length scale of the system over which the magnetic field varies strongly. 
Upon setting L ~ 10 fm and a « O.SCem^ ^ 6 MeV at T = 300 MeV, the diffusion time is 
about fo ~ 3 fm. However, as argued by Mclerran and Skokov [98], in this case the condition 
u » 1 / fc ~ 1 / Id is not satisfied, so it is more realistic to solve Eq. (19) instead its diffusion-type 
simplification. 

(2) If u 7 ^ 0 and fhe magnetic Renolds number R,„ = LUa ^ 1 (the magnetic Renolds number 
quantifies the ratio of the convection term over the "diffusion term"), we can approximately keep 
only the convection terms in Eqs. (17)-(18). This corresponds to the ideally conducting limit. The 
equations such obtained are 

() T3 

— = V X (u X B), (22) 

E = -vxB. (23) 


It is well-known that Eq. (22) leads to the/rozen-zn theorem for ideally conducting plasma, i.e., the 
magnetic lines are frozen in the plasma elements or more precisely the magnetic flux fhrough a 
closed loop defined by plasma elements keeps constant. Thus the decay of the fields are totaly due 
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to the expansion of the QGR To see the consequence of Eqs. (22)-{23), we assume for simplicity a 
initial Gaussian transverse entropy density profile 

s{x,y)=s„exp|^-^-|i^, (24) 

where ax,y are the root-mean-square widths of the transverse distribution. They are of order of 
the nuclei radii if the impact parameter is not large. For example, for Au + Au collision at RHIC, 
Ux ^ Uy 3 fm for b = 0, and Ux ^ 2 fm, fly ~ 3 fm for b = 10 fm. By assuming the Bjorken 
longitudinal expansion, 

Vz = (25) 


one can solve the ideal hydrodynamic equations for transverse expansion and obtain [99], 

Cg 

Vx = -|xf. 


'jy = -^yt, 


(26) 

(27) 


where Cg = \/dP/de is the speed of sound. Substituting the velocity fields into Eq. (22), we can 
solve out B{t) analytically. For example, the By(f) at r = 0 is give by 


By{t,0) 




(28) 


This is just manifestation of the frozen-in theorem, because the areas of the cross section of the 

2 

QGP expands according to ^ exp ^(f^ — fg) in x — z plane, thus the total flux across the x — z 
plane is a constant. Setting fl;c ~ fly ~ 3 fm and Cg ~ 1/3, we see from Eq. (28) that for f < 5 fm 
By{t) decays approximately as By{t) oc {to/t)By{to) — much slower than the 1/f^-type decay in 
the insulating case discussed in last subsection. 

So far, we discussed two special cases of Eqs. (17)-(18) which permit analytical treatments. It is 
desirable to solve the most general equations in accompanying with the hydrodynamic or kinetic- 
equation simulation for the, i.e. the fluid velocity and temperature, of the fireball. But up to today, 
this is not done yet. 


H. Electromagnetic fields in other collision systems 

In recent years, RHIG has also run heavy-ion collisions of nuclei other than gold, for example, 
the Cu + Au and U + U collisions. The EM fields in these collision systems were also studied [9, 12- 
14,100]. We here give a brief summary of these studies. 

The EM fields in U + U collisions were studied thoroughly in Ref. [9]. The uranium nu¬ 
cleus, unlike the Au or Pb nucleus, has a highly deformed prolate shape. But this shape deforma¬ 
tion does not bring much effect to the event-averaged magnetic fields. As simulated in Ref. [9], the 
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U + U collisions at = 193 GeV produce an event-averaged magnetic field just lightly smaller 
than that in Au + Au collisions at = 200 GeV, see the left panel of Fig. 8 . The readers can find 
more information, especially those related to the event-by-event fluctuation, in Ref. [9]. 




FIG. 8. The event-averaged EM fields in U + U and Cu + Au collisions. Figures are taken from Refs. [9,13]. 


The Cu + Au collisions are geometrically asymmetric: both the charge number and the total 
size of the gold nucleus are much larger than that of copper nucleus. Thus the Cu + Au collisions 
may be able to produce nonzero electric fields along the in-plane Au-to-Cu direction even after 
the event average. The numerical simulation presented in Ref. [13] found that: 

(1) The strength of the event averaged magnetic field (which is along the —y direction) in Cu + 
Au collisions at 200 GeV is comparable to that in Au + Au collisions. 

(2) There is a strong event averaged electric field pointing from the Au nucleus to Cu nucleus 
which is at the order of one tn^, see the right panel of Fig. 8 . 

(3) The azimuthal angle of the electric field, xpE, has a strong back-to-back correlation with xpi, the 
first harmonic angle of fhe participants; this is the same as that in Au + Au collisions [ 8 ]. The new 
feature is that in noncentral Cu + Au collisions there is a clear positive correlation between xpE 
and xp 2 signaling a persistent in-plane electric field. More details can be found in Ref. [13]. 

Very recently, there were interesting proposals for the novel effects of this in-plane E field in 
Cu + Au collisions, for example: the E field can lead to a directed flow vi splitting between posi¬ 
tively and negatively charged hadrons [12,14], the presence of the in-plane E field may strongly 
suppress or even reverse the sign of fhe charge-dependent correlation (see Eq. (78) for its defi¬ 
nition) [13], and Cu + Au may serve to test the chiral electric separation effect [101]; more relevant 
discussions are given in Sec. IV. 

III. ANOMALOUS TRANSPORTS IN P- AND C-ODD QUARK-GLUON PLASMA 

The strong EM fields may induce a variety of novel effects to the quark-gluon plasma, among 
which we will focus in this section on the ones that are deeply related to the topology and sym¬ 
metry of QCD and QED. It was found that, in addition to the normal electric current driven by 
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E field, there can emerge three new currents in P- and C-odd regions in QGP as responses to the 
applied EM fields. They are the chiral magnetic effect (CME) [2, 15], the chiral separation effect 
(CSE) [16,17], and the chiral electric separation effect (CESE) [18]. Thus the complete response of 
the P- and C-odd QGP to external EM field can be expressed as: 




(29) 


/ O'VV CTVA 

\aAv 

where Jy and Ja represent vector and axial currents and cr's are corresponding conductivities. 
The Ohm's law and the conductivity Uyy ^ is physically well understood, so we will not discuss 
them. In this section, we will focus on the CME, CSE, and CESE which are anomalous in the sense 
that their appearances are closely related to the topologically nontrivial vacuum structure of QCD 
and the axial anomaly. We will discuss their experimental consequences in next section. 


A. Chiral magnetic effect 


(1) What is the chiral magnetic effect? — The CME is the generation of vector current by ex¬ 
ternal magnetic field in chirality-imbalanced (P-odd) medium. Historically, the CME has been 
studied through different theoretical approaches in a number of contexts ranging from astro¬ 
physics [102], condensed matter systems [103-111], QCD physics [2, 15,112-121] , to holographic 
models [122-134]. The recent reviews of CME are [27, 29, 30, 135]. 

The CME can be neatly expressed as 


Jv = o'vaB, 

„2 


CTvA = 


2 . 71 ^ 


Pa, 


(30) 

(31) 


for each specie of massless fermions with charge e, where the current Jy is defined by ]y = 
e{xpj^‘xp) and }Ia is a parameter that characterizes the chirality imbalance of the medium. The }iA 
is commonly called aixal or chiral chemical potential although it actually does not conjugate to 
any conserved charges of the fermions; we will discuss its meaning later. Eor QGP, the total CME 
current is obtained by adding up all the light quark's contributions and the CME conductivity 
should be ayA = ^cV-aY^// (2/1^) with c\j the charge of quark of flavor / and Nc = 3 the 
number of color. 

Erom Eq. (30), we can first recognize that CME is P odd because Jy (a P-odd quantity) and 
B (a P-even quantity) transform differently under parity. Thus CME can occur only in P-odd 
medium characterized by finite }Ia- Second, the CME is C even as Jy and B are both C-odd. 
Third, the CME is T-even as both Jy and B are T-odd this is also evident from the fact that the 


^ We will specifically refer to the vector current Jy as the Lf(l) electric current J unless otherwise stated. Then uyy is 
then just the usual electric conductivity which we denoted by a in last section. 

® Here, "T" stands for "time reversal". 
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FIG. 9. The triangle diagram for CME. It is the same diagram that produces the axial anomaly. 


CME conductivity Uya, as expressed in Eq. (31), is temperature independent. The time-reversal- 
evenness of the CME conductivity indicates that the emergence the the CME current is a non- 
dissipative phenomenon [29, 121]. (Note that the usual electric conductivity ayv is T-odd and 
thus can generate entropy.) 

To get an intuitive understanding of Eq. (30), let us consider a system with unequal numbers of 
right-handed (REI) and left-handed (LH) u quarks (for example, consider that Nr > Nr with Nr/r 
the total number of RH/LH quarks) subject to uniform magnetic field. We know that the Landau 
quantization has the property that the lowest Landau level permits only one spin polarization 
which minimizes the excitation energy of quarks but is highly degenerate with a degeneracy 
factor proportional to the total magnetic flux. Thus if the magnetic field is strong enough so that 
this degeneracy factor is larger than Nr or Nr, all the u quarks are confined to the lowest Landau 
level on which their spins are totally polarized to be along the direction of the magnetic field. 
Now the RH u quarks will prefer to move along their spin direction that is the direction of the 
magnetic field; while the LEI u quarks will prefer to move opposite to their spin direction which 
is opposite to the direction of the magnetic field. Because the number of REI quarks are larger 
than the number of LH quarks, the overall effect of the motion of quarks will be to generate a net 
current of u along the magnetic field. Eor u quark, similar argument leads to a net current of u 
moving opposite to the magnetic field or a net electric current along the direction of the magnetic 
field This is the intuitive picture of the CME. Note that if the magnetic field is not strong 
enough so that the higher Landau levels are also occupied, then each higher Landau level will 
contain equal numbers of spin-up and spin-down quarks and the CME current of spin-up quarks 
will exactly cancel the CME current of spin-down quarks for each higher Landau level. Therefore 
the higher Landau levels do not contribute to the total current — only the lowest Landau level is 
responsible to the CME. 

(2) How to derive the CME? — The emergence of the CME, Eq. (30), is due to the axial anomaly 
in QED sector which couples the vector current Jy to the magnetic field B and the axial chemical 

® One should be noticed that the antipartile of a RH-chirality (RH-helicity) massless fermion is of LH-chirality (RH- 
helicity) 
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potential jiAf see Fig. 9. There are a variety of methods to derive Eq. (30) from microscopic quan¬ 
tum field theory [2, 15,102-104,115,117,119,120,136, 137], mesoscopic kinetic theory [138-152], 
to macroscopic hydrodynamic approach [22,153-158]. Here we pick up one of these derivations 
given by Fukushima, Kharzeev, and Warringa [115] because it is elementary and easy to see the 
relation between CME and the lowest Landau level and the axial anomaly. 

Let Q be the thermodynamical potential of fermions of charge e in a magnetic field B = Bz at 
finite vector and axial chemical potentials, fiy = {Br + Fl) /2 , Ba = {Br ~ /2. In the noninter¬ 

acting limit, Q can be written as 



where n runs over all the Landau levels and s is over spins, En,s is the dispersion relation of the 
fermions 


[sgn(pz) (p2 + 2 m |eB| )H2 + s^^]2 + ^2^ 


(33) 


and oi„^s is a degenerafe consfant given by a„^s = 1 — ^0n^s,-sgn{eB) which accounts the fact that 
only one spin state occupies the lowest Landau level. The current Jy can be obtained through 
differentiation of Q with respect to vector potential A^, Jy = dCi./dAz\A^=Q, which owing to gauge 
invariance is equivalent to Jy = edCl /in the integrand, 

J^ dpz dFf5,S 




2k 


s=±„=o 


[1 upi^EfiA Bv) Bv)] / 


(34) 


where A is a ultraviolet cutoff that guarantees the finiteness of the calculation at the intermediate 
steps and goes to infinity at the end of the calculation, and np(x) = (exp (v/T) -|- 1)“^ is the 
Fermi-Dirac distribution. It is then straightforward to find that 

e|eB| 


]v = 


/'lTr\2 Xj Xj [^n,s{pz — A) En^ipz — A)] 

s=±n =0 


y y a. 

{InYhh 


(A^ -6 2M|eB|)^^^ -6 s}Ia - ((A^ -h 2n|eB|)^^^ - s}Ia) 


e^BA 

2k^ 


B. 


(35) 


This is just Eq. (30). This derivation shows that: (1) The CME is due to the ultraviolet surface 
integral and is unaffected by infrared paramefers, like T, pv, etc. (2) Only the lowest Landau 
level contributes to CME, reflecting the fact that only the lowest Landau level permits a touching 
node of opposite chirality. This touching node is known to be responsible for the Ll,^(l) axial 
anomaly [103, 159]: the applied E and B fields pump the fermions at LH-chirality branch to 
RH-chirality branch at the touching node at a rate [e^/{2n^)]E- B. 

We emphasize again that although the derivation here is for non-interacting system and relies 
on Landau quantization picture, the CME conductivity Kva = ^Ba/ (2/1^) is actually fixed by the 
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axial anomaly equation and thus universal no matter how strong the interaction between fermions 
is (Recall that the axial anomaly equation itself is universal in the sense that it does not receive 
perturbative correction from scattering between fermions, a result usually referred to as Adler- 
Bardeen theorem). This is particularly supported by the derivation of CME based on holographic 
models [122-134] which intrinsically describe strongly coupled system, where although the whole 
setup is very different from the above derivation and other calculations based on perturbation 
theory, the CME conductivity is shown to be given by the same universal result. 

(3) How can QGP be chiral? — The appearance of CME requires a nonzero fiA which char¬ 
acterizes the strength of the chirality imbalance of the system. Then the question is: how can 
the QGP generate a net chirality imbalance? To answer this question, following the argument in 
Ref. [2], let us first consider the vacuum state of the SU{3) gauge theory. To make the energy min¬ 
imized, the vacuum must satisfy the condition = 0 (G";, is the field strength tensor) which 
requires the gauge field to be pure gauge: Afi{x) = ig^^U~^{x)dfiU{x) with U{x) G SU{3). 
Working in temporal gauge Ao{x) = 0 and by noting that for any time-independent gauge 
transformation, doU(x) = 0, the temporal gange fixing condition is unchanged, 0 = Aq{x) —> 
U^^{x)Ao{x)U{x) + ig^^U~^(x)doU(x) = 0, one can realize that the vacuum is described by a 
time-independent Ai{x) which is a pure gauge potential 

Ai{x) = ig~^U~^{x)diU{x). (36) 


If we impose the boundary condition U{x) —> constant at |a;| —> oo (see, for example. Refs. [160- 
163] for relevant discussions about the boundary condition), the gauge transformation U{x) de¬ 
fines a map from (R^ with the infinity identified as an ordinary point) to SU{3) which is char¬ 
acterized by a winding number riu, G 713 (561(3)) = Z, 

Hr, = J d^xe^ihr[{U-%U){U-^djU){U-^dkU)]. (37) 

This winding number is a topological invariant as can be checked by smoothly deforming U{x). 
Thus the U{x)'s corresponding to different rii^ are topologically distinct in the sense that then 
cannot be smoothly deformed into each other without passing through gauge field configurations 
whose field strengths are nonzero. In other words, the U{x)'s of different define multiple 
degenerate vacua (called the 0-vacua) separated by finite energy barriers. 

On the other hand, one can categorize all the gauge field configurations in to topologically 
distinct classes characterized by different values of the following topological invariant, 

<? = (38) 

where Ga^ = je^^^P‘^Gpg-. This q is called the (second) Chern number of configuration A and is 
always an integer [164, 165]. It is straightforward to show that under the condition G^y —)■ 0 at 
infinity (f —>• ±oo or |a;| —)■ oo). 


q = 


24 / 1 ^ 


di:peP''P‘^tT[{U~^dyU){U-^dpU){U-^dM)]r 


(39) 
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where is a surface at infinity in four-dimensional spacetime and U{x) G SU{3) satisfies 
U{t,x) —> U±{x) for t —> ±oo and U{t,x) —)■ constant for \x\ —)■ oo [160-163]. Thus, after several 
steps of manipulations, 

q = = co) - n^{t =-co). (40) 


This means that the gauge field configuration which goes to pure gauge at infinity and has finite 
q can induce a transition from vacuum of winding number = — oo) to another vacuum of 
winding number = oo). At zero temperature, such gauge field configurations are called 
instantons [166] and they are responsible for the quantum tunneling through the energy barrier 
between vacua of different winding numbers [167-170]. 

The high energy barrier (~ Aqcd ~ 200 MeV) between two vacua suppresses the instanton 
transition rate exponentially, but at high-enough temperature, the transition between different 
vacua can also be induced by another, classical, thermal excitation called sphaleron [171, 172] 
which, instead of tunneling through the barrier, can take the vacuum over the barrier. In elec- 
troweak theory, sphaleron transitions cause baryon number violation and may be important for 
the cosmological baryogenesis [173,174]. In QCD, the existence of the sphaleron configurations at 
finite temperature enormously increases the transition rate [175-184]. At very high temperature, 
the perturbative calculation of the sphaleron transition rate gives [2,176-184] 


Tsph ~ i‘^sNcfT\ (41) 

while the strong coupled holographic approach gives an even larger rate Tgph = {g^Nc)^/{256n'^)T^ [185]. 
Thus at high temperature, the sphaleron transition rate can be very large. This provides a machin¬ 
ery of generating P and CP odd bubbles in QGP (note that a transition process from a topologically 
trivial vacuum to a topologically nontrivial vacuum violates P and CP symmetry as is evident 
from the integrand of q). 

Now if we integrate the axial anomaly equation in the QCD sector (for massless quarks) 


where (/ is over all massless flavors) is the axial current, we see that 


(42) 


Na(f = co) - NA{t = -oo) = 2q, 


(43) 


where Na = f (PxJ^{x) is the total chirality or axial charge. This demonstrates that a topolog¬ 
ically nontrivial gauge field configuration can create or annihilate the total chirality of fermions, 
and thus if the QGP contains a (sufficiently large) domain in which q is finite we would expect 
that it finally will contain unequal numbers of RH and LH quarks or anti-quarks even if initially 
NA{i = — oo) = 0. This is how QGP can become chiral. We note here that the probabilities of 
generating positive chirality and negative chirality are equal which means that over many col¬ 
liding events in heavy-ion collisions the averaged chirality should vanish. What remains after 
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event average is the chirality fluctuation rather the chirality itself and any measurement of the 
chirality-imbalance effects should be on the event-by-event basis, see Sec. IV. 

(4) What is axial chemical potential? — There is a conceptual problem in interpreting as 
the axial chemical potential for fermions: as when we talk about the chemical potential we always 
need to associate to it a conserved quantity but we know that the axial charge J of fermions 

is in general not conserved when the fermions are coupled to gauge fields. Furthermore, if fiA is 
really a chemical potential conjugate to a conserved axial charge of fermions, then Eq. (30) would 
imply a persistent electric current even at global equilibrium. But the appearance of CME current 
in equilibrium violates the gauge invariance in the QED sector as discussed in Ref. [29], see also 
Refs. [109, 124,126, 135] for relevant discussions. Thus the physical meaning of }Ia is actually quite 
confusing, and there have been a number of relevant discussions in literature [27, 29,109,124,126, 
135, 186, 187]. According to these studies, the CME should be considered as a non-equilibrium 
phenomenon which vanishes when the system is in global equilibrium and, correspondingly, }Ia 
shouldn't be regarded as a fixed chemical potential of the fermions in equilibrium. It is more 
appropriate to view }Ia ^s the rate of the time changing of fhe 6 field (see below), }Ia = h is 
the 6 parameter that describes the state of fhe system. 

To reveal the relation between }Ia and the 9 field lef us consider the 0-vacua of QCD. The effect 
of the 0-vacua can be encoded into a 0-term in QCD Lagrangian, 

^QCD = -^G^^Gr + E^/[^y'(a, - igA,)]xPf - 0|^G;,Gr, (44) 

where is the gluon field. For massless quarks the 0-term does not give observable consequence 
because it can be rotated away by 11 , 4 ( 1 ) transformation, but if we promofe the 0 angle to be 
spacetime dependent (a pseudoscalar 'axion' background), then it will show important physical 
consequences. Perform the path integral over the quarks: 

Z[A,e] = j [dxp][dTp]e'^‘^-^\ (45) 

The 0-term in the action can be eliminated by performing a local 11,4(1) transformation to quark 
fields and by using the Fujikawa method: 

xPf{x) ^ 

xpf{x) i/i/(x)yG)75. (46) 

The resulting Lagrangian reads 

Gqcd = - jG^yr + Ey - igA, + i(9,0)75)]i/^/. (47) 

^ / 

From this Lagrangian we can identify }Ia = dtO, that is, jia is the time derivative of the 0 field. 
Now let us consider the QED sector. Then the following Lagrangian should be added to Gqcd^ 

CQED = -lF^vF^^-AXr 


(48) 
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where is the photon field and Jy — with qjr the charge of quark of flavor f. The 

above Ua(l) transformation will then induce a new term in the the final effective Lagrangian 


^QCDxQED = “ -G"^Gr + - igA^ + 1(3^9) J 5 )]lpf - A^ify + -OFf'^'P^y, 

(49) 


where 

''=N«E^' (50) 

The EM sector of this Lagrangian is the Maxwell-Chern-Simons or axion Lagrangian. The equa¬ 
tions of motion derived from this Lagrangian is called the axion electrodynamics [188], which 
read (see e.g. Ref. [189]) 


V-E = J^ + kV 9 ■ B, 
dE 

V X B —— = Jy — k{}IaB + V9 X E), 


V-B = 0, 

dB 

V X E -h -7:r- = 0- 


(51) 

(52) 

(53) 

(54) 


Lrom the second equation, we can identify that k}IaB is the CME current which plays the same 
role as the applied current Jy. 

We see that the out-of-equilibrium nature the CME is already indicated by the identification 
= dt9- if we allow the 6 angle to be the correct physical quantity specifying the physical state of 
the system, then a constant f /,4 would already indicate a linearly growing 6 in time. Then what is 
the counterpart of CME in equilibrium? This can be guessed by integrating both sides of Eq. (30) 
over time and we can write the result down as (for constant magnetic field): 

P = 9—^B A constant, (55) 

where P is the electric polarization vector. This equation describes a magnetic-field induced 
charge polarization through the axion-EM-field coupling, which can appear in insulating matter 
even in equilibrium and is called the topological magentoelectric effect [190,191]. In particular, re¬ 
cently found time-reversal-invariant topological insulators fulfill the topological magnetoelectric 
effect with a special value of 6, i.e., 6 = n. 


B. Chiral separation effect 

(1) What is the chiral separation effect? — The chiral separation effect (CSE) is the dual effect 
to the CME in the sense that the former is expressed from the latter by interchanging the vector 
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and axial quantities [16, 17, 192], In formula, the CSE is given by 



Ja = o-aaB, 


(56) 

(57) 


where the axial current is defined by = e{xpj^'y 5 ip), and }iv is the vector chemical potential 
corresponding to the global IJv(l) symmetry (can be considered as baryon chemical potential) 


Thus the CSE represents the generation of the axial current along with the external magnetic field 
in the presence of finite vector charge density parametrized by the vector chemical potential f/y. 
Unlike the CME, the CSE is P-even and T-even but C-odd, the appearance of CSE does not require 
a parity-violating environment. Similar with the CSE, the T-even nature of CSE signals that it is 
not dissipative. 





FIG. 10. The Feynman diagram for CSE. 


The intuitive picture of the emergence of the CSE at finite vector density is the following. Imag¬ 
ine a fermionic system with more charged massless fermions than anti-fermions. At extremely 
large magnetic field pointing to up direction, supposing the fermions are positively charged, then 
all the fermions are confined in the lowest Landau levels and their spins are aligned upward. Then 
the up-moving fermions carry positive helicity while the down-moving fermions carry negative 
helicity; because the helicity of a massless fermion is equal to its chirality, we see the fermions 
contribute an axial current along the direction of the magnetic field. Similarly, it is easy to see that 
the anti-fermions contribute an axial current opposite to the direction of the magnetic field. Be¬ 
cause the system contains more fermions than anti-fermions, a net axial current moving upward 
is generated. 

(2) How to derive the CSE? — Similar with the CME, the generation of CSE is also a conse¬ 
quence of the axial anomaly, see Eig. 10. Comparing Fig. 10 to Fig. 9 one can again recognize the 
dual relation between CSE and CME. The expression (56) has been already indicated in the axial 

^ It is not necessary to define and jiy this way. We can also define the axial current as (like in 

last subsection), i.e., without the e factor, and }iy as the electric chemical potential. Or we can also define the axial 
current as and jiy as the baryon chemical potential, but then ctaa will read aAA = ejiy/Irp- 
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anomaly equation. 


dtfA + V-Ja = ■ B. (58) 

Assuming steady state and noticing that we can write eE = V}iv, this suggests that Ja = 
[e^/for uniform B. A more rigorous derivation along this line can be found in, for 
example. Ref. [17] in which the surface integral are more appropriately handled. Other deriva¬ 
tions and various aspects of the CSE are studied in Refs. [16, 17, 111, 192-200]. We here just stress 
one fact that like the CME, only the lowest Landau level is responsible to the arising of CSE. 
Intuitively, this is because the higher Landau levels are occupied by equal numbers of spin-up 
and spin-down fermions (supposing B is along the up direction), it is easy to see that they give 
opposite contributions to and thus cancel out (Similar argument applies to anti-fermions as 
well). 

(3) The chiral Magnetic Wave. — The duality between CME and CSE allows the existence 
of interesting collective modes called chiral magnetic waves (CMWs). This can be seen in the 
following way. We can express the CME (30) and CSE (56) in terms of RH and LH currents and 
corresponding chemical potentials: 


e 

Jr = 

where 

1 

Jr = ^{Jv + Ja)> 
1 

Jl = 2^Jv ~ Ja)/ 

Hr = Hv + Ha/ 

Hl = Hv- Ha- 


(59) 

(60) 


(61) 


In terms of RH and LH quantities, the continuation equations read. 


aJs/L + V ■ Jr/l = 0, 


(62) 


where and are the densities of RH and LH fermions. Now consider a small departure from 
equilibrium so that we can write and }Ir^l as }Ir^l = Ho A ^Hr,L/ Jrl — ^o + ^Jr i where 
5j\i^,5}iR^i are all small numbers. Substituting Eqs. (59)-(60) into Eq. (62) and keeping linear 
terms in 3}i and 3J^), one can obtain 


dtSjl + 





B ■ VSJl = 0, 




B ■ W5]l = 0, 


(63) 


(64) 
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where Xr = ^Jr/^P^r and Xl = dJi/djii are susceptibilities for RH and LH chiralities. These are 
two wave equations expressing collective modes arising from the coupled evolution of CME and 
CSE. They are just the CMWs [19] (see also Ref. [201] for a derivation of CMW based on the kinetic 
theory), one propagating along the direction of B and another opposite to B with wave velocities 

Vr = 


4^71^ Xr' 


Vl = 


(65) 


4ti^Xl ' 

Some comments are in order. (1) If in addition to B there is also a finite E such that E ■ B 0, then 
the continuation equations are anomalous and the right-hand side of Eqs. (62) and thus the wave 
equations (63) and (64) should be added terms ±[e^/ {27t^)]E ■ B. These terms pump chirality 
into the system and provide sources to the wave equations. The CMW velocities are not affected 
by the anomalous terms. (2) In general, the microscopic scattering among constitute particles will 
lead to diffusion of the RH and LH charges. In this case, we should add the diffusion terms, 
— DV/^/^, where D are diffusion constants, into Eqs. (59) and (60). The presence of diffusion 
renders the dispersion relation of CMW to receive a —iDk^ term with k the wave number of the 
CMW, and thus dampens the mode of wave number \k\ ^Vr/r/D. 


C. Chiral electric separation effect 

(1) What is the chiral electric separation effect? — The CME and CSE are both induced by 
magnetic field. As we have seen in Sec. II that heavy-ion collisions can generate strong electric 
fields as well. Thus one may wonder whether the electric field can induce an anomalous transport 
phenomenon. Such an anomalous transport indeed exists, as first found in Ref. [18] in weakly 
coupled QED plasma and later in weakly coupled QCD plasma [202], and is called the CESE. It is 
also confirmed in some strongly-coupled holographic setups [203, 204]. 

The CESE predicts the generation of axial current in external electric field when the system has 
both finite vector and axial densities characterized by corresponding chemical potentials jiy and 
}iA respectively, 

Ja = (T'avE, ( 66 ) 

^AV = XeWpiA, (67) 

where Xe is not universal and depends on the microscopic scattering processes. Unlike the CME 
and CSE, the generation of CESE is not directly related to triangle anomaly in the QED sector. 
Actually, the one-loop diagram for CESE has four legs which is not responsible for the triangle 
anomaly, see Eig. 11. To understand why }ivHA, let us examine the CPT properties of 

Ja and E: is P-even, C-even, and T-odd while E is P-odd, C-odd, and T-even. Thus aAv 

must be P-odd, C-odd and T-odd. This means that aAv should be proportional to odd powers of 
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fly and when jiy/T and JIa/T are small (as indeed in high-energy heavy-ion collisions) the 
leading contribution would be proportional to }Ia}I-v/T'^- On one hand, the T-odd property of cfav 
suggests that the CESE may be dissipative, but on the other hand Cav is not positively definite 
because the sign of jiy (and }Ia) can be eifher positive or negative suggesting CESE to be a non- 
dissipative phenomenon. Detailed study [202] revealed that the dissipation property of CESE is 
very subtle: CESE can affect the entropy production and in this sense it is dissipative in nature, 
but the CESE always appears in accompany with the usual electric conduction, it alone does not 
generate entropy. This is why aAv can be either positive or negative. 


wvw 

liy 

E 

wvw 



Ea 

ww 


FIG. 11. The Fe 5 mman diagram for CESE. 


(2) The CESE conductivity. — The CESE conductivity aAv depends on the microscopic interac¬ 
tion. It can be calculated along the same line as the computation of fhe usual electric conductivity 
(Tw In Ref. [18] the CESE conductivity for hot QED plasma was calculated by using the Kubo for¬ 
mula within the Hard Thermal Loop perturbation theory [84, 205-208]. At the leading-log order 
in the coupling constant e the result reads 


QED 


20.499 


EuEa T 

e2ln(l/e)' 


( 68 ) 


where one should note that the axial current Ja in the present paper is e times the axial current 
defined in Ref. [18] so that there is factor e difference between (Tav here and that in Ref. [18]. The 
similar perturbative calculation was extended to hot QCD plasma later, and the result for the 
two-flavor case that u and d quarks are massless reads [202] 


QCD,2f 

^AV 


14.5163 Tr/(Q,Q,4) 


}iv}iA e^T 
T2 g^\n{l/g)' 


(69) 


Note that in this case the axial current is defined as = e{tpjf‘'Y 5 QAip) and Je — e{ip'Y^Qerp) with 
Qe = diag(2/3, —1/3) and Qa a generic matrix in flavor space; while the chemical potentials 
}iv and }iA are corresponding to Lly(l) and Ua{1) currents. Eor the case with u,d, and s quarks 
massless, the result is [202] 

= (70) 
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The CESE conductivity was also studied in certain holographic model [203] and it was found 
that cx; which was shown to be surprisingly robust even for large 

jiv/T and jiA/T. 

(3) The collective modes due to CESE. — With the CME, CSE, and CESE found, we can write 
down the complete response of the chiral medium to external EM field as (here we focus on QED 
plasma for simplicity) 

Jy = CTyvE + (TyaB = dyvE + (75}IaE, (71) 

Ja = ctavE + ctaaB = XeHvHAE + dsjivB, (72) 

where U 5 = e^/ (27r^) and Xe is the numerical factor in the CESE conductivity. Because dvA ^ }Iai 
(Taa ^ }>-Vr and (Tav HvI^a, we see that in the presence of external EM field, the vector and axial 
densities mutually induce each other and get entangled in a nontrivial way. As a consequence, the 
coupled evolution of small fluctuations of the two densities can induce collective modes in a way 
similar to that the energy density fluctuation in usual hydrodynamics lead to sound wave. The 
first such example is the CMW which we have discussed in last subsection. The CESE introduces 
nonlinearity (in the }iv}i-A term) in to the equations and make the problem more complex [18]. 

To see how the collective modes arise, we consider static and homogeneous external E, B 
fields, and study the coupled evolution of the small fluctuations in vector and axial charge den¬ 
sities. Eirst, we can write the total EM field as £?tot = E + 5E and Stot = B -\- 5B with 3E cx. 
and SB oc Jy are induced fields due to the vector density and current fluctuations. Second, tak¬ 
ing into account that the chemical potentials are small, we can replace the chemical potentials 
by the corresponding charge densities: jiy^A = ^v,aJv a' where the OLy^A are the inverse suscep¬ 
tibilities defined as ocy^A = ^V-v,aI^1% a- Third, we expand dyy as dyy = Uq -|- cr 2 {}iy + }i\) and 
omit higher order terms in the chemical potentials. Then by linearizing the continuity equations 
/v A + ^ ■ Jv,A = 0, one can obtain 

dfSjl + dQSfy + (Ts0La{B ■ V)Sj\ 

+2a24ny{E ■ V)Sfy + 2a2a^AnA{E ■ V)5i\ = 0, 
dt5i\ + asOLv{B ■ V)5]^ + Xe‘>^v<^Any{E ■ 'V)df^ 

+Xe‘>cvCiAnA{E ■ V)Sjy = 0, (73) 

where ny and tia are small uniform background vector and axial densities and Sjy ^ ^ nv,A are 
fluctuations on top of them. 

Without loss of generality, we suppose B is along z-axis, i.e., B = Bz while E = Ee. The 
dispersion relation obtained from Eq. (73) is 

[m - Vy{e-k)] ± ^\J[m - V-{e-k)f + ^A^{k), (74) 

where v± = Va with = 2 a 20 LynyE and Va = Xe‘^v‘^AnyE, and 

A^{k) = [cr 5 a.AB{z ■ k) + 2a20c\nAE{e ■ k)] x [cr 5 «.yB{z ■ k) + XeOi-vOiAnA^i^ ■ k)]. (75) 
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To reveal the physical content of the above dispersion relation, we consider three special cases. 

(1) B = Bz and E = 0. 

Then Eq. (74) reduces to a; = ±^^{v^kz)^ — {cro/2)^ — i{crQ/2) with speed = a^-^fcLy^B. When 
^ uo/2 it represents two propagating modes (jO ~ ±v^kz — f(uo/2): these are generalized 
CMWs which reduce to the CMWs discussed in last subsection when ctq = 0 and ocy = a a- When 
— ^^o/2 the two modes are damped. 

(2) E = Ez, B = 0, and ny = 0. 

In this case, we find two modes from (74) 

CV = ±\J{Vekz)^ - ((7o/ 2)2 - f( ctq/2) (76) 

with Vt; = aAnA\/'^o' 2 Xe<x-vix-A^- Similar to the CMWs, when v^k^ S> (Tq/2 there are two well- 
defined modes co ~ Evgkz — i(ao/2) from CESE that propagate along E field and can be called 
the Chiral Electric Waves (CEWs). They become damped when Vekz < o'q/2. 

(2) E = Ez, B = {), and nA = 0. 

In this case, the vector and axial modes become decoupled, and Eq. (74) leads to 

C0y{k) = Vykz - iao , 

CVAik) = Vakz. (77) 

They represent a "vector density wave" (VDW) with speed Vy — 2a20c\nyE that transports vector 
charges along E field but will be damped on a time scale ~ 1 /no and a propagating "axial density 
wave" (ADW) along E with speed Va = Xe‘>‘-v‘>‘-AnyE without damping. 

These collective modes may be able to transport chirality and charges in heavy-ion collisions 
and can lead to specific charge azimuthal distributions from which the CESE may be detected 
experimentally. We now turn to discuss the experimental status of the searches of the anomalous 
transports. 


IV. EXPERIMENTAL SEARCHES OF THE ANOMALOUS TRANSPORTS 

QCD dynamics does not violate parity globally which means that in heavy-ion collisions after 
average over many events there will not remain a finite axial chemical potential }Ia — The Ha 
just fluctuates from event to event. As a result, the P-odd anomalous transport phenomena can 
only be measured on the event-by-event basis. We now discuss the event-by-event observables 
for CME, CMW, and CESE. 

A. The charge-dependent azimuthal correlation and chiral magnetic effect 

The CME transports charges along the direction of B. In noncentral heavy-ion collisions, the B 
direction is, on average, perpendicular to the reaction plane, thus the CME is expected to induce 
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a charge separation with respect to the reaction plane. After the hydrodynamic evolution of the 
fireball, this out-of-plane charge dipolar configuration in coordinate space can be converted to 
special distribution of the final charged hadrons in momentum space which can be detected by 
utilizing specifically designed charge-dependent azimuthal correlations. One such correlation has 
been proposed by Voloshin [209] (see also Ref. [112]), 

7a/? = {cos{(pa + (pfi- 2Yrp)), (78) 

where the indices a. and /3 denote the charge of the hadrons, (pa and (p^ are the azimuthal angles 
of hadrons of charge a and j 6 respectively, Yrp is the reaction plane angle, and (■ ■ •) denotes 
average over events; See Fig. 12 for illustration. It is easy to understand that a charge separation 

with respect to the reaction plane will give a positive 7^ _and a negative 7++ or 7 _In real 

experiments, the measurements are done with three-particle correlations where the third hadron 
(of arbitrary charge) is used to reconstruct the reaction plane, 

{C0s{(pa + (pji-2(p3)) 5^27a/?/ (79) 

where V 2 is the elliptic flow. 



FIG. 12. Illustration of a t 5 rpical noncentral collision. Modified from Ref. [210]. 


Experimental results for 7 ^^. — The first measurement of 7 ^^ was officially reported by STAR 
Collaboration at RHIC for Au + Au and Cu + Cu collisions at = 200 GeV in 2009 [211,212], see 
Fig. 13 (PHENIX Collaboration at RFIIC also reported equivalent measurement [213, 214]). From 
several aspects, these experimental results are consistent with the expectation of the CME: (1) 
Away from the central collision region, the data shows a clear positive opposite-sign correlation 
and negative same-sign correlation as the CME predicted. The magnitude of the measurement 
is in agreement with the theoretical estimate by Kharzeev et al. [2, 112]. (2) The observed cor¬ 
relation increases from zero at the most central collisions to the peripheral collisions, consistent 
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with the fact that the magnetic field increases as the centrality grows. (3) The date shows that the 
opposite-sign correlation is smaller that the same-sign correlation at a fixed centrality, in agree¬ 
ment with the fact that the opposite-sign particles are emitted in opposite directions and thus 
their correlation suffers from quenching because at least one of the two correlated particles needs 
to propagate through the opaque medium. Recently, the STAR Collaboration officially reported 
the measurement of at different beam energies [210, 215]. The main message of the new mea¬ 
surement is that the CME-induced signal is reduced at lower energy; in particular, for energy 
smaller than 20 GeV, the same-sign and opposite-sign correlations become to overlap with each 
other. This is consistent with the fact that at low temperature the QCD topological transition rate 
is strongly suppressed and thus the CME is attenuated; when the temperature is so low that the 
chiral symmetry is spontaneously broken, the large mass effect will strongly diminish the axial 
chemical potential arid thus the CME. Note that the magnitude of the magnetic field is roughly 
proportional to \/s, thus one may think that at lower energy the magnetic field is less efficient to 
transport charges. Elowever, this is not likely the case, as the lifetime (as defined in Sec. II E) of 
the magnetic field due to the spectators is inversely proportional to \/s which compensates the 
reduction of its magnitude when \/s decreases. Recently, the AEICE collaboration at LHC pub¬ 
lished their measurement of for Pb + Pb collisions at = 2.76 TeV [216]. The result shows 
that the magnitude of the observed correlation is very close to that at RHIC in spite that the 
magnetic field at LHC is about much larger than at RHIC. 



FIG. 13. The measured charge-dependent correlation y^^ by STAR Collaboration. This figure is from 
Ref. [211]. 


Let us briefly summarize the current experimental results for the charge-dependent correlation 
that are consistent with the CME interpretation. 

• In high-energy collisions (RHIC @ = 200 GeV and LHC @ y^ = 2.76 TeV): 

(1) Positive opposite-sign correlation and negative same-sign correlation for nonzero centrality 
are clearly seen. 

(2) The observed magnitude of the charge-dependent correlation is at the same order as CME 
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predication. 

(3) The correlation grows as the centrality increases. 

(4) The magnitude of the opposite-sign correlation is relatively smaller that that of the same-sign 
correlation at the same centrality. 

(5) The magnitude of the correlation at given centrality is nearly independent of the collision 
energy. 

• In STAR beam-energy scan, \/s = 7.7 — 200 GeV: 

(6) The difference between same-sign and opposite-sign correlations decreases when collision 
energy drops and finally nearly disappears for energy less than ~ 20 GeV. 

The background effects. — In the STAR publications [210-212,217], they also reported results 
for another, reaction-plane independent, charge-dependent correlation, 

3afi = {cos{(pcc-4>fi)). (80) 

If the CME is the only source for the charge-dependent correlations, then one would expect that 
= 3 _>0 while <i-|_< 0 which is easy to understand from the out-of-plane charge separa¬ 
tion picture. However, the data shows the opposite: ^3 _<0 while > 0. This tells us 

that the measured charge-dependent correlations contain contributions from background effects 
other than CME and these background contributions are likely at the same order of magnitude 
of (or even larger than) the CME contribution [218]. A number of studies have been conducted 
aiming to identify and quantify possible background effects [13, 116, 219-234], but so far there 
still lack a quantitative and reliable way to estimate and subtract those backgrounds. Detailed 
discussions on possible background effects can be found in the review papers [235, 236]. We here 
will just briefly discuss two possible background effects, namely, the transverse momentum con¬ 
servation [222, 223, 227] and the local charge conservation [222, 225, 226]; they are considered to 
comprise the most important backgrounds. 
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FIG. 14. Illustration for the transverse momentum conservation and its contribution to 7 ^^. 


The momentum conservation enforces the sum of the transverse momenta of all the parti- 
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cles to be zero at any time during the evolution of the system (here "transverse momentum" 
means the momentum component transverse to the beam direction). Thus if we observe a parti¬ 
cle of transverse momentum p^, the transverse momentum conservation (TMC) requires all the 
other particle must have a total transverse momentum —; see Fig. 14 (left-panel) for an illustra¬ 
tion. This builds an intrinsic back-to-back two-particle correlation (here "intrinsic" means that this 
correlation is irreducible in the sense that it carmot be described by single-particle distribution) 
and it is evidently charge-blind. Suppose on average each particle has a transverse momentum 
of the same magnitude and the total multiplicity is N ^ 1. Then the TMC requires that 
cos{(j)i — Yrp) = sin{(j)i — Yrp) = 0 for each event where (pi is the azimuthal angle of 

particle i. Then we can obtain that 


^ TMC _ „ TMC _ / COs(^, -I- (pj 2 Yrp) ^ 

Tss - Tos - \--/ 

^ ^ [E; cos{4)i - Yrp)]^ - [£,■ sin(^,- - Yrp)]^ - cos(2^,- - 2 Yrp) ^ 

1 v? 

= -j;^^(cos(2cp, - 2 Yrp)} - 


irTMC _ j;TMC _ fj) ^ 

%s — ^OS — \ / 

^ ^ [Li cos{(pi - Yrp)]^ -F [Zi sinj^i - Yrp)]^ - ^ 


L//; 


1 


N-1 


1 

'n' 


(81) 


(82) 


where SS and OS denote "same-sign" and "opposite-sign" and we have used the assumption 
that N 3> 1. Although the TMC itself is reaction-plane independent, the elliptic flow V 2 knows 
the reaction plane; furthermore the elliptic flow also "measures" the centrality. The TMC always 
gives negative contributions to both 7 and and it becomes more important for more peripheral 
collisions because fhe multiplicity decreases and Vz increases as the centrality decreases. 



FIG. 15. Illustration for the local charge conservation and its contribution to 7 ^^. 


The TMC can give a large amount of contribution to 7 ^^ and as shown by AMPT simu¬ 
lation [230, 237] (This can also be roughly estimated by choosing vz ~ 0.1 and N ~ 1000 for 
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moderate centrality bins at RHIC which gives ~ 10^^ and ^™c ^ being at the same 
order of magnitude as the STAR results [211, 212]). In order to evidence the contributions of 
CME, it is useful to subtract the TMC contribution from 7 . This can be implemented by taking the 
difference between opposite-sign and same-sign correlations, 

^7 = ^( 7 +-- 7 ++- 7 —)• ( 83 ) 


Although this correlation does not receive contribution from TMC, it gets contribution from the 
local charge conservation (LCC) which is charge dependent. The schematic cartoon for LCC is 
depicted in Fig. 15. Suppose that in the coordinate space the charge is conserved locally in the 
sense that a cell of certain size ~ R contains zero total charge. If we pick up one such cell which 
moves along certain azimuthal direction, then upon freeze-out, the hadrons emitted from this cell 
will contain equal numbers of positive and negative charges and they collectively move along 
the same direction as the original cell. In such a way, the LCC induces a near-side opposite- 
sign correlation but no same-sign correlation. Let us now make a rough estimate of the LCC 
contribution in an idealized setup. Suppose at the freeze-out the system is consist of Nc non¬ 
overlapping cells with each containing M/2 positive and negative hadrons, thus N = MNc with 
N the total charged hadron multiplicity. Let the cell c move along the azimuthal angle (pc and 
assume all hadrons inside c emit along the same direction cpc. We can then express the opposite- 
sign correlations as (the same-sign correlations due to LCC are zero) 


T re 

7+- 


^LCC 


Lt,qi = + Lj,qj=- COs{(pi + (pj - 2 Yrp) ^ 

^ N+N_ ^ 

/ Lc Liec,qi^+ Ljec,qj=- cos{(pi + (pj - 2 Yrp) ^ 

' N+N_ ^ 

M^LcCos( 2 ,^,- 2 Yrp) 

' 4N+N_ ' 

M 

Li,q,=+Lj,qj=-COs{(Pi-(Pj)^ 

^ N+N- ^ 

'Yhc YLiec,qi^+ Ylijec,qj=- , 

' n+nI ' 

M 

N' 


(84) 


(85) 


where N+ and N_ are the multiplicities of positively and negatively charged hadrons and we 
have used N+ ~ N_ = N/2 3> 1 and Nc S> 1. More thorough analysis can be found in Refs. [222, 
225, 226] in which the effects of finite size of the cell and the finite emission angle (the hadrons 
may emit from the cell c within a small relative angle {(pc — A(p/2,(pc + A(p/2)) are more accurately 
addressed by using the charge balance function. The LCC contribution is inversely proportional 
to the total multiplicity and proportional to the volume of the neutralized cell. Comparing to the 







35 


TMC we see that LCC is very strong for both 7 and 5 as it is proportional to M which is larger 
than 1 . 

Clearly, to fully understand the experimental data and to determine the CME-mduced signals, 
one has to build a realistic modeling of the charge-dependent correlations in which all the strong 
backgrounds and the CME should be taken into account on the same footing. By assuming that 
the dominant backgrounds are those driven by elliptic flow, like TMC and LCC, it is plausible to 
make a two-component modeling of the two-particle correlations as [9, 235], 

lafi = KV2Fafi - (86) 

(87) 

where the F and H represent the strengths of background effects and CME, respectively. The 
constant k is an "uncertainty factor" accounting for some practical issues such as finite kinematic 
selections in real experiments that may drive k to deviate from unity [210]. By adopting these 
relations, one may use the data for Sap, and V 2 to extract the strength factors F and H as 
functions of the centrality. Recently, the STAR collaboration reported the measurement of AH = 
Hss — Hqs [ 210 ] and found it is finite, increasing with the centrality, and insensitive to collision 
energy: features consistent with the expectation of CME. 

The U + U and Cu + Au collisions. — It was proposed that the U + U collisions can serve 
as a potential mean to distinguish the flow-driven effects from the CME in 7 ^^ [238], and such 
collisions were recently done at RHIC with some preliminary results reported in [215]. The main 
point is that, unlike the gold or lead nuclei, the uranium nucleus is highly deformed from spheri¬ 
cal shape and even in the most central U + U collisions there can be sizable elliptic flow while the 
magnetic field is expected to be very tiny [9]. Therefore, if 7 ^^ is dominated by CME, the signal 
would disappear in the most central U + U collisions; while if is dominated by flow-driven 
backgrounds like TMC and LCC, there will still be sizable signal even for very central collisions. In 
Ref. [9], by substituting the F and H factors extracted from the Au - 1 - Au data, the computed mag¬ 
netic field, as well as the V 2 data for U + U collisions [215] to Eq. ( 86 ), the authors are able to make 
prediction for A 7 in U + U collisions. Although there is still quantitative discrepancy between the 
prediction and the preliminary experimental data [215], the two show consistent magnitudes and 
trends. This suggests that the two-component scenario may be a good approximation toward a 
more reliable modeling of the charge dependent correlations. But more experimental releases and 
more quantitative simulations are definitely needed, see for example Refs. [100, 239] for recent 
attempts. 

In Ref. [13], a new possible way by using the Cu - 1 - Au collisions to test the CME was proposed. 
The main idea is the following (see Eig. 16 for an demonstration). As we have seen in Sec. IIH 
that Cu + Au collisions can generate a in-plane electric field as well as the out-of-plane magnetic 
field. This in-plane electric field can drive a in-plane charge dipole in addition to the out-of-plane 
charge dipole due to CME. As numerically shown in Ref. [13], the appearance of the in-plane 
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dipole makes both 73 $ and 70 s descend or get reversed if the electric field is much stronger than 
the magnetic field. However, this picture is based on the assumption that 70 -^ is dominated by 
EM-field induced effect; if it is dominated by U 2 -driven effects rather than EM-field effect, we do 
not expect that 7 ^^ changes too much from Au + Au collisions to Cu + Au collisions. (A plausible 
guess would be that A 7 as a function of centrality in Cu + Au collisions lie between that in Cu + 
Cu and Au + Au collisions.) Thus, measuring 7 ^:^ in Cu + Au collisions may give useful insight 
into the different mechanisms underlying 7 ^/ 3 . 







centrality 


FIG. 16. Illustration of the charge dipole induced by CME and electric field in Cu + Au collisions. (Left): 
No in-plane electric field. The dipole is driven by CME to be perpendicular to the reaction plane. (Middle): 
Weak electric field. The in-plane electric field drives the dipole to skew from the magnetic field direction. 
Both 7 os and 753 should be reduced. (Right): Strong electric field. The in-plane electric field drives a large 
in-plane component of the charge dipole, and thus the signs of 70 s and 733 will be reversed. This figure is 
from Ref. [13] 


B. The charge-dependent elliptic flow and chiral magnetic wave 

As we have seen from Sec. IIIB, the CMWs are able to transport both chirality and electic 
charge. Phenomenologically, they can induce an electric quadrupole in the QGP [78, 240]. This 
can be schematically seen from Pig. 17. In the heavy-ion collision, the overlapping region contains 
a small amount of vector density (inheriting from the the colliding nuclei) so that fly > 0 
Therefore the B field can induce an axial current along it via the CSE (56) and leads to a chirality 
separation with respect to the reaction plane with > 0 in one tip of the overlapping region 
and f/A < 0 in another tip. The CME in turn induces two vector currents which transport positive 
charges toward the two tips and negative charges toward the equator of the overlapping region 
and eventually form an electric quadrupole. Qualitatively, the strength of the electric quadrupole 
is proportional to }iv{eB)^. 

® Taking into account the event-by-event fluctuations, there can be events with overall piy < 0 in the overlapping 
region. The induced electric quadruple will then flip. 
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Then the next question is: what is the observable signal for this electric quadrupole? It was 
proposed in Ref. [78] that such a signal could be the splitting between the elliptic flows of 7r+ 
and 71^; see Refs. [241-244] for more advanced numerical implementation of this proposal. As 
illustrated in Fig. 17, owing to the strong in-plane pressure gradient, the electric quadrupole in 
coordinate space can be converted into a specific pattern in the momentum-space distribution 
of the final-state hadrons after freeze-out: more negatively charged hadrons move along the in¬ 
plane direction while more positively charged hadrons move along the out-of-plane direction. As 
a result, the elliptic flow of will be larger than that of tt^. This is true for events with jiy > 0; 
For a event with y < 0, the result will be that ^ 2 ( 71 +) > V 2 {tc^)- To quantify this V 2 splitting, one 
can consider the single-particle azimuthal distribution fq {(p), = ±, per rapidity. The appearance 
of the electric quadrupole modifies the net charge distribution to the following form: 

N-u - N- 

Pi^P) =/+(</’)-/-(</’) =- 2 ^^-[1-F2z;2Cos(2</)-2Yrp) -rcos(2</)-2YRp)], (88) 

where r oc (eB)^ is a parameter quantifying the strength of the electric quadrupole at given — 

N- CX lly, 

2 

^ = - #cos(2(/>-2Yrp) [p(</))- p(<^)|b=o], (89) 

V 2 is the elliptic flow when the electric quadrupole is absent, and N+, N_ are the multiplicities 
of positively and negatively charged hadrons. Suppose the quadrupole is small so that the total 
distribution is unchanged, 

f+i4>)+f~W = + ^^2 cos (2./) - 2 Yrp)] . (90) 

From Eq. (88) and Eq. (90), we can obtain the single particle distribution of charge q as 

fq{<p) ~ ^ [1 + 2^2 cos(2</) - 2Yrp) - qAaj{2(p - 2Yrp)] , (91) 

where A^h — (N+ — N_)/ (Ny + N_) is the net charge asymmetry parameter and we have kept 
only terms up to linear order in A^r. Erom Eq. (91), we immediately see that the elliptic flow 
becomes charge dependent: 


V2{n^) =V2T 




2 


(92) 
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We give two comments here. (1) Obviously, the presence of the electric quadrupole in the QGP 
should affect not only the elliptic flow of charged pions but also the elliptic flows of other charged 
hadrons. The reason why we finally take into account only pions, according to the argument in 
Ref. [78], is because other abundant charged hadrons, the proton, antiproton, kaons, etc, suffer 
from large difference in the absorption cross section in hadronic matter at finite baryon density 
which may strongly mask or even reverse the the effect of the electric quadrupole. (2) We should 
also notice that in Eq. (88) we assumed that the quadrupole moment is perpendicular to the reac¬ 
tion plane. This is not exactly true because, as we saw in Sec. HE, B orientation fluctuates event- 
by-event; by taking into account this azimuthal fluctuation of the B, the parameter r should be 
proportional to ((eS)^ cos[2(i/jb — Yrp)]) rather than {(eB)^) [8]. 

Recently, the V 2 splitting of charged pions was measured by STAR Collaboration [215, 245, 246] 
and the data is in agreement with the CMW predictions in several aspects (See Eig. 18 for the ex¬ 
perimental result for the centrality dependence of the slope parameter r at = 200 GeV): 

(1) The measured V 2 difference ^ 2 ( 71 ^) — ^ 2 ( 71 ^) is linear in with a positive slope as predicted 
by CMW; see Eig. 18. 

(2) The magnitude of the measured slope r versus centrality can be fitted by the CMW estimation 
by adopting the magnetic field simulated from the event-by-event calculations with field lifetime 
Tb ~ 4 — 5 fm. 

(3) The data for the slope parameter r as a function of centrality shows a maximum in midcentral¬ 
ity or midperipheral collisions which is also consistent with the CMW calculation [246]. 

(4) The slope parameter r displays no obvious trend of the beam energy dependence for 10 — 60% 
centrality at ^/s = 20 — 200 GeV. This is consistent with the fact that the ability of fhe magnetic- 
field in transporting charges in a tiny rapidity interval is approximately energy independent be¬ 
cause the strength of the magnetic field due to spectators is proportional to v^s while the its life¬ 
time is proportional to 1 / 

(5) At = 11.5 and 5.5 GeV the data shows that the slop parameters are consistent with zero — 
a feature in agreement with CMW picture because the arising of CMW requires chiral symmetry 
restoration which is not expected to be the case at very low ^/s. 

Although the experimental data shows very encouraging features in line with the CMW ex¬ 
pectations, one cannot exclude the possibility that other effects may contribute to the observed V 2 
splitting pattern as well. Indeed, the data (see the left panel of Eig. 18) shows clearly a positive 
^2 () — ^2 () at Ajh = 0 which is not understood within (but also not in contradiction with) 
the CMW picture. A variety of possible background effects were proposed, for instances, the 
out-of-plane electric-field can induce out-of-plane electric currents which can produce a charge 
quadrupole similar with that produced by CMW (Note that this electric-field induced effect exists 
even at A^h = 0) [7, 79], the electromagnetic chirality E ■ B distributes in the transverse plane like 
a chiral dipole which when combined with triangle anomaly and CME can also induce a charge 
quadrupole [79], the mean-field potential from the hadronic as well as partonic phases acts differ- 
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FIG. 18. (Left) The vi splitting between n~ and 7r“*" as a function of charge asymmetry. (Right) The mea¬ 
surement of the slop parameter r by STAR Collaboration as well as the UrQMD and CMW predictions. 
Figures are from Ref. [246]. 


ently on particles and anti-particles and produces a difference between and 7r+ [247, 248], 
the effect of baryon stopping, when assuming that the Vi of transported quarks is larger than that 
of produced ones, can induce a larger V 2 of than 7r+ [249], the effect due to the rapidity cut 
in the data analysis when combined with the local charge conservation in rapidity can lead to 
splitting of charged pions [250], the isospin effect combined with viscous hydrodynamics may 
also show a Vi splitting [251]. It should be noted that none of the above backgrounds can success¬ 
fully explain all the features of the experimental data. It is quite plausible that different effects 
are entangled to produce the measured pattern. To distinguish different effects and separate the 
CMW contribution, we certainly need more detailed measurements such as the measurement of 
the V 2 . splitting versus transverse momentum and rapidity, the measurement for other hadrons 
like kaons, and the measurement of splitting in other harmonic flows like v-^. 


C. The possible observables for chiral electric separation effect 

Because the CESE is driven by electric field, in order to observe the CESE in heavy-ion colli¬ 
sions one must have a collision system where the electric field is generated. As seen in Sec. IIB, 
owing to the event-by-event fluctuations of the proton locations inside the nuclei, heavy-ion col¬ 
lisions do generate strong electric field even in the center of the collision region. However, as we 
discussed in Sec. IIE, so-generated electric field is not correlated to the matter geometry (the re¬ 
action plane or participant plane), and thus invisible in any reaction plane dependent observable. 
To have an electric field that correlated to the matter geometry, we need to consider asymmetric 
collisions. The Cu -i- Au collision is such a collision where a persistent electric field in generated 
with the direction from Au nucleus to Cu nucleus. Thus let us consider the possibility of detection 
of the CESE in the Cu -i- Au collisions [18]. 

In Cu -I- Au collisions, supposing that the created hot QGP contains both vector and axial 
charge densities from fluctuations and topological transitions, then the E field will lead to both 
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FIG. 19. Illustration of the charge configuration due to the in-plane charge separation effect due to usual 
electric conduction, CESE, and CME in Cu + Au collisions. The figure is from Ref. [18]. 


an in-plane charge separation via usual Ohm's conduction and an in-plane chirality separation 
via CESE. The resulting in-plane axial dipole will then further separate charges via CME along 
the magnetic field which is in the out-of-plane direction, and cause an approximate quadrupole at 
certain angle ipc in between in- and out-of-plane, see Fig. 19 for demonstration. After the collective 
motion and freeze out of the system, this spacial distribution of charges can cause corresponding 
single-particle distribution of charged hadrons in momentum space, which can be expressed as 

fqi(p) = ^{1 + 2yicos(^ - Yi) -F2u2Cos[2(^ - Y 2 )] +2qdECOs{(p - ips) 

+2xqdB cos{(p — xps) + '2-Xd^B cos[2((^ — ipc)] + higher harmonics}, (93) 

where q = dz is the charge of the particle, Nq is the number of hadrons with charge q per rapidity, 
2 are the usual directed and elliptic flows at .B = E = 0, Yi 2 are corresponding harmonic 
angles, tpc is the CESE angle (as shown in Fig. 19, positive/negative charges fly out at the angle 
±ipc and Expc + tc) which lies in (0, 7i/2), dE and dg characterize the strength of the dipole induced 
by E and B, respectively, and ;\; = ± is a random variable accounting for the fact that in a given 
event there may be a sphaleron or an anti-sphaleron transition resulting in > 0 or < 0 via 
triangle anomaly. The event average of x should be zero. Equation (93) represents four kinds of 
flows: the normal directed flow and elliptic flow (the second and third terms), and the normal 
electric conduction along tpE (the fourth term), the pure CME along tps (the fifth term), and the 
combined CESE and CME (the sixth term, which we will simply call CESE term). 

From the single-particle distribution (93), we first observe that the directed flow of charge q 
acquires a contribution from the electric field: 

vi{q) = vi +q{dECOs{xpE-'^ 1 )) (94) 

where (■ • •) denotes event average. Note that the CME term does not contribute to Ui {q) because 
the event average of x is zero. Thus in Cu + Au collisions, the measured directed flow is likely to 
be charge dependent. This may provide us a possibility to extract the electric conductivity of the 
produced matter by measuring the difference between the directed flows of positive and negative 
hadrons because dg is proportional to the electric conductivity. This idea has been exploited in 
Refs. [12,14]. 
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In Ref. [101], two possible observables for CESE were propoesed. The first one is the charge- 
dependent two-particle correlation 

= (cos[2((/)fl: + 2Yrp)]), (95) 

where cpa stands for the azimuthal angle of hadrons of charge a with a = ± and the average is 
taken over events. We will set Yrp = 0 for simplicity. It is easy to find from Eq. (93) and the con¬ 
figuration shown in Fig. 19 that the same-sign correlation ^gs ~ (cos(4t/?c)) should be in general 
different from the opposite-sign correlation ^os ~ 1 while the purely dipolar charge distribution 
(like the in-plane charge separation effect (in-plane CSE) due to usual conduction by E or CME, 
the fourth and fifth terms in Eq. (93)) does not contribute to Thus the correlation (95) carries 
information of the quadrupolar distribution due to CESE. However, in real experiments, may 
receive contributions from background effects related to the elliptic flow V 2 because, if we turn 
off all the anomalous effects, ~ ^ 2 - Another background may be the transverse momentum 
conservation [222, 223, 227] which causes a contribution cx v^/N with 174 the fourth harmonic 
flow and N the multiplicity. Thus, just like what has been done for CME observeble 7 ^^, we can 
take the difference = ^os ~ Css to subtract these charge-blind backgrounds. But there may 
remain another backgrounds, for example, due to the local charge conservation [222, 225, 226] 
(which leads to Cos oc v^/N and ^gg ~ 0 and thus AC oc V 4 /N) or due to the chiral magnetic 
wave induced quadrupole. These backgrounds again need to be more carefully explored and 
subtracted in any real experimental measurements. The AMPT simulation for AC is shown in 
Fig. 5 of Ref. [101] for Cu + Au collisions from which one can clearly see that when the CESE is 
present the increasing trend of AC from central to peripheral collision are present. 

The second observable for CESE is the charge-dependent event-plane angles. Define AY = 
(IY^ — Y 2 I), where Y^ and Y 2 ^ are the second harmonic angles of the event planes reconstructed 
by final positively charged hadrons and negatively charged hadrons. Then from Fig. 19 one can 
find that AY ~ 2 (t/ 7 c). The AMPT simulation for AY is given in Fig. 6 of Ref. [101]: As expected, 
once the CESE effect happens, a visible AY is present which shows a linear dependence on the 
centrality while if for the initial settings that do not include CESE, the nonzero AY is absent for 
the whole centrality window, making AY a good observable for CESE. 


V. DISCUSSIONS 

In this review, we have discussed the recent progresses in understanding the properties of the 
strong electromagnetic (EM) fields generated in heavy-ion collisions and the anomalous transport 
phenomena induced by these EM fields. The discussions have shown that heavy-ion collisions, 
the only mean on earth to create the quark-gluon plasma, can also generate extremely strong EM 
fields. These EM fields provide us the possible monitor of the non-trivial topological structure of 
the quantum chromodynamics (QCD) via the anomalous transport phenomena. We have focused 
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on several types of the anomalous transports, namely, the chiral magnetic effect (CME), chiral 
separation effect (CSE), chiral electric separation effect (CESE), and the corresponding collective 
modes, e.g., the chiral magnetic wave (CMW), emerging from the mutual induction of the vector 
and axial vector charges and currents through them. 

We have also discussed the ideas and the current status of the experimental searches of these 
anomalous transports. Hitherto, the experimental measurements show very encouraging fea¬ 
tures that are qualitatively consistent with the expectations of the CME and CMW, although many 
background effects may give significant contributions to the observables as well, making the inter¬ 
pretation of the experimental measurements very challenging. More efforts are certainly needed 
before the final conclusions can be made. On the theoretical side it is very desirable to make real¬ 
istic studies of the CME and other anomalous transports in heavy-ion collisions based on kinetic 
simulations or (magneto)hydrodynamics by consistently encoding the time evolution of the EM 
fields, the event-by-event generation of the axial chemical potential and its spatial distribution 
and time evolution, etc. 

We keep all the discussions as pedagogical as possible so that this review can be accessed by 
readers who are not familiar with the topics that we have covered. We give literature at proper 
places where deeper and more technical treatments can be found. This review will be poten¬ 
tially useful for graduate students who want to begin their study in the area of EM-field induced 
anomalous transport phenomena in heavy-ion collisions. 

There are a number of interesting topics that are untouched so far but are closely related to 
the subjects covered by this article. Here we briefly address three of them which have attracted 
remarkable attention recently: 

(1) The noncentral heavy-ion collisions can generate fluid vorticity in the interaction re¬ 
gion [252-258]. The fluid vorticity can also induce anomalous transports, generally called the 
chiral vortical effect (CVE) [20-22] (earlier suggestion is given in Ref. [113]), which are the vorti¬ 
cal analogues of the chiral magnetic and separation effects. The corresponding collective modes 
owing to the CVE are the chiral vortical waves (CVW) [23] just like that the chiral magnetic waves 
are owing to CME and CSE. Phenomenologically, the CVE can induce baryon number separation 
with respect to the reaction plane [259] which can be tested by using the baryon number depen¬ 
dent correlation = {cos{(pa + (pb — 2Yrp)) with here (pot the azimuthal angle of hadrons of 
baryon number a (=baryon or antibaryon) and the CVW can induce a V 2 splitting between A and 
A in a way very similar with how the CMW induce V 2 splitting between and 7r+. Recently, the 
STAR Collaboration reported the first measurement of which shows evident splitting between 
the same-sign (i.e., baryon-baryon or antibaryon-antibaryon) correlation and the opposite-sign 
(i.e., baryon-antibaryon) correlation for large centrality[260]. This is consistent with the expecta¬ 
tion of CVE. However, just like the case of the CME detections, the experimental signal is possibly 
masked by contributions from other effects like the transverse momentum conservation and the 
local baryon number conservation, and more efforts are certainly needed to understand the data. 
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(2) Quite recently, the triangle anomaly was reported to be realized in special condensed mat¬ 
ter systems called Weyl semimetals [261-265]. These materials are (3+l)-dimensional analogues 
of (2+l)-dimensional graphene. Their band structure permits the existence of the band-touching 
points at which the dispersion relation of the quasiparticles is approximately linear and thus ef¬ 
fectively "relativistic". Near a band-touching point the equation of motion of the quasiparticle 
can be described by Weyl equation and hence the name Weyl semimetal; the corresponding band¬ 
touching point is called Weyl point. The Nielsen-Ninomiya theorem [266, 267] requires that the 
Weyl points in the Brillouin zone of a Weyl semimetal must come in pairs with opposite chiral¬ 
ity. If two Weyl points with opposite chiralities happen to overlap with each other, the effective 
equation of motion becomes Dirac type and such materials are called Dirac semimetals [268-272]. 
The Dirac semimetal can be effectively transmuted into a Weyl semimetal by applying parallel 
electric and magnetic fields. Actually, in such a way, the authors of Ref. [273] reported the first 
observation of the chiral magnetic effect in Dirac semimetal ZrTe^- The discovery of the Weyl and 
Dirac semimetals opens a new era for the studying of the topological matter and we can expect 
that other anomalous transport phenomena, like CMW and CESE, could also be realized in Weyl 
and Dirac semimetals in future. 

(3) In addition to the Weyl and Dirac semimetals, the spin-orbit coupled atomic gases may 

provide another possibility to realize the anomalous transports in condensed matter systems. The 
synthetic spin-orbit coupling (SOC) was first successfully generated for Bose gas of afoms 
in which two hyperfine states (referred to as (pseudo)spin-up and -down states) of the atoms 
are coupled to their orbital motion in a manner described by Elamiltonian AH ~ — the 

so-called Rashba-Dresselhaus SOC Elamiltonian— where ax is a Pauli matrix, ky is y component 
of the momentum, and A is the SOC strength [274]. The same type of SOC was also realized in 
Eermi gases of ^Li [275] and [276], and the pure Rashba SOC was generated in Eermi gas 
of [277]. Very promisingly, the Weyl SOC which can be described by the Hamiltonian H = 

/ {2m) — Xcr ■ k with m the mass of atoms could also be realized [278-280]. In Ref. [Ill], based 
on the chiral kinetic theory, the author showed that the triangle anomaly can be realized in Weyl 
spin-orbit coupled atomic gases in a harmonic trap under rotation. In this case, the trapping 
potential provide an effective electric field and the rotation provides an effective magnetic field. 
As a consequence of this triangle anomaly, the currents of opposite helicities are generated in 
parallel or anti-parallel to the rotation axis which mimic the chiral magnetic effect and chiral 
separation effect. The potential experimental signature is that the chiral magnetic wave in this 
case can induce a mass quadrupole in the atomic cloud which may be detected by, e.g., light 
absorption images. 
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